For  Reference 


NOT  TO  BE  TAKEN  FROM  THIS  ROOM 


The  University  of  Alherta 
Printing  Department 
Edmonton,  Alberta 


UNIVERSITY  OF  ALBERTA 
LIBRARY 

Regulations  Regarding  Theses  and  Dissertations 


Typescript  copies  of  theses  and  dissertations  for  Master’s  and  Doctor’s 
degrees  deposited  in  the  University  of  Alberta  Library,  as  the  official  Copy  of 
the  Faculty  of  Graduate  Studies,  may  be  consulted  in  the  Reference  Reading  Room 
only. 


A  second  copy  is  on  deposit  in  the  Department  under  whose  supervision  the 
work  was  done.  Some  Departments  are  willing  to  loan  their  copy  to  libraries, 
through  the  inter-library  loan  service  of  the  University  of  Alberta  Library. 

These  theses  and  dissertations  are  to  be  used  only  with  due  regard  to  the 
rights  of  the  author.  Written  permission  of  the  author  and  of  the  Department 
must  be  obtained  through  the  University  of  Alberta  Library  when  extended  passages 
are  copied.  When  permission  has  been  granted,  acknowledgement  must  appear  in  the 
published  work. 

This  thesis  or  dissertation  has  been  used  in  accordance  with  the  above 
regulations  by  the  persons  listed  below.  The  borrowing  'library  is  obligated  to 
secure  the  signature  of  each  user. 


Please  sign  below: 


Date 


Signature 


Institution 


THE  UNIVERSITY  OF  ALBERTA 


OPTIMAL  CONTROL  OF  HIGH  SPEED  TRAINS 


by 


Lloyd  Edwin  Peppard 


A  THESIS 

SUBMITTED  TO  THE  FACULTY  OF  GRADUATE  STUDIES 
IN  PARTIAL  FULFILMENT  OF  THE  REQUIREMENTS  FOR  THE  DEGREE 

OF  MASTER  OF  SCIENCE 


DEPARTMENT  OF  ELECTRICAL  ENGINEERING 


EDMONTON,  ALBERTA 


SPRING,  1969 


Digitized  by  the  Internet  Archive 
in  2019  with  funding  from 
University  of  Alberta  Libraries 


https://archive.org/details/Peppard1969 


UNIVERSITY  OF  ALBERTA 


FACULTY  OF  GRADUATE  STUDIES 


The  undersigned  certify  that  they  have  read,  and 
recommend  to  the  Faculty  of  Graduate  Studies  for  acceptance, 
a  thesis  entitled  Optimal  Control  of  High  Speed  Trains  sub¬ 
mitted  by  Lloyd  E.  Peppard  in  partial  fulfilment  of  the 
requirements  for  the  degree  of  Master  of  Science. 


ABSTRACT 


In  recent  years  the  problem  of  providing  high  speed 
ground  transportation  for  large  urban  populations  has  received  consid¬ 
erable  attention.  In  this  thesis,  a  review  of  the  basic  concepts  of 
high  speed  train  operation  is  presented  leading  up  to  the  development 
of  a  mathematical  model  for  a  system  of  trains.  The  model  is  based  on 
a  previous  derivation  using  Newton's  Second  Law  but  is  modified  to 
include  jerk  (rate  of  change  of  acceleration)  as  a  variable. 

A  train  regulator  is  derived  from  known  theory  of  opti¬ 
mal  control  (the  linear  state  regulator  problem  with  infinite  terminal 
time)  and  the  resultant  control  system  simulated  on  an  analogue  com¬ 
puter.  System  behaviour  is  compared  with  that  of  the  previously  devel¬ 
oped  system  using  several  performance  indices  in  the  form  of  cost 
functionals . 

The  optimal  regulator  developed  in  this  thesis  appears 
to  better  meet  the  objective  that  excessive  jerk  be  minimized. 
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CHAPTER  1 


INTRODUCTION 


1.1  The  Balanced  Transit  System 

For  the  past  several  years  the  problem  of  urban  trans¬ 
portation  has  generally  been  attacked  by  the  construction  of  vast  free¬ 
way  systems  to  facilitate  the  movement  of  private  automobile  traffic. 
The  automobile  continues  to  be  a  desired  form  of  transportation  because 
of  its  comfort,  convenience  and  privacy.  Most  cities  also  have  a  well 
developed  bus  system  which  in  general  utilizes  the  same  roads  and  free¬ 
ways  as  private  and  commercial  automobile  traffic. 

[1] 

As  shown  by  Michaels  ,  modern  high-speed  rail  rapid 
transit  is  the  most  efficient  and  economical  method  known  for  moving 
large  numbers  of  people  quickly  and  comfortably.  Further,  a  "balanced" 
transit  system  with  rail  transportation  as  its  "core"  and  including  bus 
and  automobile  traffic  can  be  the  most  advantageous  scheme  for  urban 
centres.  This  scheme  would  result  in  less  congestion  and  less  air  pol¬ 
lution  and  hence  benefit  user  and  non-user  alike. 

An  interesting  comparison  can  be  made  between  the  effi¬ 
ciency  of  various  forms  of  transportation  in  the  following  manner.  We 
will  assume  that  on  a  lane  of  freeway  traffic  there  is  an  average  of 
1.5  persons  per  car  and  2.5  sec.  spacing  between  cars  (180  feet  at 
50  mph) .  Then  the  lane  is  capable  of  moving  2160  passengers  per  hour. 
Now  let  us  assume  a  train  of  5  cars,  each  capable  of  carrying  100  pas¬ 
sengers,  A  train  separation  of  1.5  min.  would  result  in  a  flow  rate  of 
20,000  people  per  hour  or  the  equivalent  of  about  9  lanes  of  freeway. 
Further,  in  the  case  of  buses,  for  36  second  spacing  and  50  passengers 
per  bus  the  flow  rate  would  be  5,000  people  per  hour.  Although  a 


_ 


2 


number  of  assumptions  are  made  in  arriving  at  the  above  flow  rates,  they 
do  indicate  that  a  higher  efficiency  is  attainable  by  a  train  system. 

In  high  density  areas  train  operation  could  be  made  even 
more  efficient  by  operating  with  smaller  headways  (the  distance  from 
the  rear  of  one  train  to  the  front  of  the  following  train  is  termed  the 
headway) .  In  the  above  discussion  spacings  were  given  in  units  of  time 
which  allows  calculation  of  headway  in  units  of  distance  for  any  veloc¬ 
ity,  A  further  discussion  of  headways  is  given  in  Chapter  2. 

From  the  above  it  is  clear  that  high  speed  rail  lines 
provide  transportation  of  highest  efficiency  in  areas  where  the  density 
of  passengers  is  high  enough  to  warrant  short  headways.  Most  of  the 
existing  high  speed  rail  systems  were  built  before  the  use  of  the  auto¬ 
mobile  had  become  widespread.  Hence  they  were  restricted  to  urban  areas 
af  high  population  density  such  as  New  York  City.  Today  in  areas  such 
as  Los  Angeles ,  the  population  density  is  not  as  high  as  in  New  York 
but  the  use  of  the  automobile  and  bus  as  feeder  systems  for  a  centrally 
located  rail  system  can  contribute  to  an  equally  efficient  operation. 

Such  a  balanced  system  could  find  applications  in  most  North  American 
cities  with  populations  exceeding  400,000,  This  would  include  such 
Canadian  cities  as  Winnipeg,  Vancouver  and  Edmonton. 

The  last  consideration  which  we  will  mention  here  is 
that  of  cost  of  the  system.  The  cost  of  adding  rail  lines  to  existing 
freeway  systems  need  not  be  high.  For  example,  in  Chicago,  a  rail  line 
was  constructed  in  the  median  of  the  4  lane  Eisenhower  expressway  at  an 
increase  in  total  expenditure  of  only  10%.  In  many  areas,  sections  of 
existing  rail  line  could  be  converted  for  use  in  a  balanced  transit 
system  thus  reducing  the  cost  considerably.  Certainly  in  many  areas, 
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congestion  has  become  so  serious  a  problem  as  to  make  cost  of  a  superior 
system  less  of  a  consideration. 

One  primary  requirement  for  the  successful  operation  of 
the  high  speed  systems  discussed  above  is  the  development  of  safe,  reli¬ 
able  and  effective  automatic  control  systems  for  controlling  vehicle 
speed,  acceleration  and  spacing  while  running  at  the  most  efficient  head¬ 
way.  Passenger  comfort  becomes  a  prime  consideration  and  results  in 
several  constraints  being  placed  on  any  envisioned  control  system.  These 
will  also  be  considered  in  Chapter  2. 

1.2  Two  Operational  Systems 

In  this  section  two  of  the  several  existing  high  speed 
rail  systems  will  be  discussed  briefly  so  as  to  indicate  the  specifi¬ 
cations  of  these  systems  in  regard  to  top  speeds,  headways,  etc. 

The  Bay  Area  Rapid  Transit  (BART)  System  of  San  Francisco 
is  perhaps  the  best  example  of  a  modern  high  speed  rail  system  as 
envisioned  in  the  previous  section.  Scheduled  to  begin  revenue  opera¬ 
tions  in  1970,  the  system  has  the  following  design  standards: 

1.  The  trains  run  on  separated  rights  of  way  including  aerial 
structures,  subways,  and  underwater  tubes  which  free  the 
entire  system  from  grade  crossings. 

2.  Propulsion  is  provided  by  a  1000-volt  d-c  vehicle  propulsion 
and  power  supply  system  based  on  a  1000-volt  d-c  contact 
rail  and  solid-state,  chopper-controlled  traction  motors. 

3.  Straight  line  performance  specifications  are  acceleration 
from  zero  to  50  mph  in  20  seconds  in  contrast  to  the  usual 
train  standard  of  20  seconds  from  zero  to  25  mph.  The  top 


cruising  speed  is  80  mph. 


, 


-  A  - 

A.  A  headway  of  about  90  seconds  is  maintained  by  the  automatic 
train  control  system  which  is  also  in  charge  of  all  train 
movements  such  as  acceleration,  speed  control  and  deceleration. 
Provision  is  also  made  for  manual  operation. 

The  primary  principle  of  regulation  used  in  the  BART  scheme  is  that  all 
trains  moving  over  a  given  section  of  track  will  be  moving  at  the  same 
speed.  This  will  be  a  prime  requirement  of  the  control  system  des¬ 
cribed  later  in  this  thesis. 

Another  train  system  which  has  been  in  operation  for  a 
considerably  longer  time  than  the  BART  System  is  the  New  Tokaido  Line 
of  the  Japanese  National  Railways.  Commercial  passenger  service  was 
opened  between  Tokyo  and  Osaka  on  October  1,  196 A  and  has  since  oper¬ 
ated  successfully  with  a  traffic  load  of  up  to  150,000  passengers  per 
day. 

Since  it  is  an  inter-city  system,  the  specifications  of 
the  Tokaido  System  differ  considerably  from  those  of  the  intra-city 
BART  System  described  earlier. 

1.  The  maximum  operating  speed  is  130  mph  with  a  minimum 
headway  of  5  min. 

2.  Deceleration  requires  5,600  feet  to  slow  from  130  mph  to 
100  mph  (first  stage)  and  5,250  feet  to  slow  from  100  mph 
to  zero  (second  stage).  Including  pre-braking  coasting, 
it  takes  11,600  feet  to  stop  from  a  speed  of  130  mph. 

The  resultant  maximum  deceleration  rate  is  1.6  mph/sec. 

3.  The  control  system  operates  as  follows:  Each  section  of  the 
line  has  a  specified  speed  depending  on  the  nearness  of 
stations,  curves  and  the  presence  of  other  trains  on  the 
line.  On  entering  a  section  at  other  than  the  required 


- 


5 


speed,  acceleration  or  braking  is  applied  automatically „ 

The  normal  length  of  a  controlled  block  of  track  is  10,000 
feet . 

4,  In  the  future,  operational  speeds  of  up  to  160  mph  are 

[2] 

envisioned  which  according  to  Fuji!  is  close  to  the 
maximum  practical  speed  limit  for  a  conventional  railway 
of  standard  gauge. 

Due  to  this  practical  speed  limit  much  attention  has  been 
paid  in  the  literature  to  new  and  often  radical  methods  of  propulsion 
and  guidance.  As  we  are  primarily  interested  in  the  control  of  con¬ 
ventional  systems,  these  will  not  be  described  here, 

1,3  Scope  of  This  Thesis 

In  the  following  chapters  a  high  speed  train  regulator 
will  be  developed  using  existing  theory  of  optimal  control  (the  state 
regulator  problem),  A  mathematical  model  of  the  system  is  developed 
from  Newton's  Second  Law.  The  work  reported  is  basically  an  extension 
of  the  work  of  Levine  and  Athans^^^  and  differs  in  the  following 
respects.  Jerk  (time  rate  of  change  of  acceleration)  is  taken  into 
account  by  differentiating  the  equation  of  motion  for  each  train. 

Several  performance  criteria  (in  the  form  of  cost  functionals)  are 
used  to  evaluate  the  optimal  control  system's  performance,  A  compari- 
tive  study  of  the  system  of  Levine  and  Athans  and  that  developed  in 
this  thesis  is  made  by  simulating  both  control  systems  on  an  analogue 


computer. 


_ 
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CHAPTER  2  PROBLEM  FORMULATION 

2.1  Description  of  the  System  and  Objectives 

The  system  of  trains  to  be  considered  is  represented  in 
figure  2-1.  There  are  N  trains  in  the  system  of  which  only  a  typical 
section  of  three  is  shown  in  figure  2-1.  The  mass  of  the  1th  vehicle 
is  m^  and  its  scheduled  position  at  time  t  is  given  by  z^Ct).  The  head¬ 
way  between  vehicles  is  h  and  the  length  of  each  vehicle  is  L. 

In  general,  there  are  two  ways  of  regulating  the  motion 
of  a  string  of  vehicles.  Either  each  vehicle  can  be  made  to  regulate 
its  velocity  and  position  with  respect  to  an  internal  on-board  reference, 
or  a  communications  system  can  provide  each  vehicle  with  a  knowledge 
of  the  velocity  and  position  of  every  other  vehicle  in  the  string  thereby 
enabling  the  regulation  of  differences  in  velocity  and  position  of  adja¬ 
cent  vehicles.  It  is  fairly  obvious  that  the  first  type  of  regulation 
would  not  prevent  collisions  in  the  case  of  uncontrolled  motion  of  one 
or  more  vehicles  in  the  string  (as  in  the  case  of  a  de-railment) .  The 
second  type  of  regulation  will  control  the  headway  and  for  this  reason 
we  will  deal  primarily  with  regulation  of  this  type, 

2.2  Assumptions  and  Constraints 

The  following  assumptions  are  made  in  formulating  the 
mathematical  model  for  the  system: 

1.  It  will  be  assumed  that  the  system  of  vehicles  is  running 
on  a  flat  section  of  straight  guideway.  This  assumption 
is  important  in  the  derivation  of  the  equation  of  motion 
for  each  vehicle. 

2.  The  operation  is  assumed  to  take  place  far  enough  between 
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Figure  2-1  THE  PHYSICAL  SYSTEM 
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stations  that  controlled  braking  or  acceleration  is  not 
required. 

3.  The  vehicles  making  up  the  train  are  assumed  to  be  of 
equal  mass.  This  will  in  practise  be  realized  only  when 
each  vehicle  is  carrying  an  equal  number  of  passengers. 
Normally,  there  will  not  be  too  great  a  deviation  from 
this  condition  when  the  system  is  operating  in  an  efficient 
manner . 

We  are  now  in  a  position  to  consider  the  constraints 
imposed  on  the  system  which  must  be  taken  into  account  when  formulating 
the  problem. 

1.  Velocity  -  The  operating  speed  of  the  system  of  trains  is 
self-constrained  by  the  limitations  of  the  physical  equipment 
such  as  propulsion  and  guidance  devices  and  thus  will  not 

be  considered  in  the  formulation  of  the  problem. 

2.  Acceleration  and  Deceleration  -  In  the  interest  of  passenger 
comfort  we  must  constrain  the  acceleration  forces  generated 
by  the  controller  to  be  less  than  some  maximum  figure. 
Friedlander ^ ^  has  suggested  a  naximum  limit  equal  to  1.4 
times  the  pull  of  gravity.  However,  it  is  felt  that  under 
emergency  conditions  (i.e.  to  avert  a  collision)  that  higher 
forces  could  be  used. 

It  is  interesting  to  note  that  the  maximum  acceleration- 
deceleration  rates  of  the  BART  and  Tokaido  Systems  are  both 
less  than  0.1  g.  which  is  well  below  that  experienced  in  the 
private  automobile  (0.8  g.  for  a  good  disc  brake  system). 


- 
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3.  Jerk  -  Time  rate  of  change  of  acceleration  force  or  jerk 
force  is  more  objectionable  to  passengers  than  the  accel¬ 
eration  force  itself. 


FIGURE  2-2  A  POSSIBLE  ACCELERATION  TRAJECTORY 


Figure  2-2  shows  a  possible  plot  of  acceleration  vs  time  for 

correction  in  velocity  of  a  train  in  the  system.  The  time 

delay  At  in  application  of  acceleration  a^  is  due  to  the  sum 

of  all  time  delays  in  the  system  including  inertia  of  the 

vehicle.  Depending  on  At,  the  maximum  jerk  given  by 

(da/dt)  can  be  excessive.  In  the  private  automobile,  the 

max  • 

driver  is  in  control  of  the  applied  jerk  at  all  times  and 
generally  applies  acceleration  forces  (except  in  emergency 
cases)  in  such  a  manner  as  to  make  the  maximum  jerk  accept¬ 
able  to  him  yet  accomplishing  his  required  velocity  change 
in  a  reasonable  length  of  time.  This  is  exactly  what  we  would 
like  to  accomplish  with  our  train  regulator. 


10 


4.  Headway  -  For  a  number  of  reasons  this  parameter  must 

impose  a  rigid  constraint  on  our  system.  As  indicated  in 
figure  2-1,  headway  is  defined  as  the  distance  from  the 
rear  of  one  train  to  the  front  of  the  following  train.  If 
collisions  are  to  be  avoided  the  minimum  headway,  h^ ,  is 
chosen  so  that  under  emergency  conditions  a  vehicle  is  able 
to  come  to  a  stop  without  colliding  with  the  preceding 
vehicle.  The  maximum  line  capacity  is  given  by 


C 


v 


h„  +  L 


0 


(2.1) 


where  v  is  the  cruising  velocity  and  L  the  vehicle  length, 
c 

r  6 1 

Hojdu  et  alL  give  the  minimum  stopping  distance 
especially  in  the  case  of  regenerative  braking  as 


v 


- 


_  f 


gdg  +  Atv 


0  0  a(a) 


(2.2) 


where  v^  is  the  velocity  before  deceleration  and  a(a)  expresses 
deceleration  as  a  function  of  the  train  velocity  (a  is 
a  dummy  variable).  The  total  communication  and  reaction 
time  delay  is  given  by  At. 

Figure  2-3a  shows  an  idealized  trapezoidal  decel¬ 
eration  trajectory.  The  corresponding  velocity  trajectory 
is  given  in  figure  2- 3b.  By  expressing  the  deceleration  in 
terms  of  the  velocity  at  any  time  it  has  been  shown ^ 
by  applying  2.2,  the  minimum  stopping  distance  is 
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X 
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c  +  JL  ^  v  4-  Atv 
2a  2  da  c  c 

-1  dt 


(2.3) 


FIGURE  2-3 


DECELERATION  AND  VELOCITY  BRAKING  TRAJECTORIES 
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If  a  fictitious  constant  emergency  deceleration  rate 
a0<aM  ^e^^-nec^  then  we  can  write 


x 


0 


+  Atv 

c 


(2.4) 


where  a^  reflects  the  effects  of  maximum  allowable  jerk. 
If  we  neglect  the  time  delay  At  we  can  write  the  minimum 
headway  as 


,  k>0 


(2.5) 


In  the  case  of  a  train  de-railment  or  track  blockage 
and  with  k<l  we  can  expect  to  have  collisions.  We  shall  term 
this  type  of  emergency  condition  a  "catastrophic  failure". 

For  the  string  of  vehicles  separated  by  the  minimum  headway 
h^ ,  the  Nth  vehicle  will  be  able  to  stop  without  colliding 
with  the  preceding  car  in  the  case  of  a  catastrophic  failure 
if 

NhQ  +  (N-l)L  >  xQ  +  (N-l)L  (2.6) 

as  shown  in  figure  2-4. 
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FIGURE  2-4 


A  CATASTROPHIC  FAILURE 
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Equation  2.6  can  be  reduced  to 


(2.7) 


Then  the  number  of  trains  that  will  collide  in  the  case 
of  a  catastrophic  failure  is 


n  =  N-l  =1-1  ,  0<k<l 

c  k 


(2.8) 


2.3  Mathematical  Model  of  the  System 

Given  the  assumptions  of  section  2.2,  Levine  and  Athans 
describe  the  motion  of  the  ith  train  in  a  string  of  N  trains  by  the 


[3] [4] 


second-order  differential  equation 


Vi  +  8±(V  =  fi(t) 


(2.9) 


obtained  by  applying  Newton's  Second  Law.  The  mass  of  the  ith  train 
is  given  by  m_^  and  its  displacement  by  x_^.  The  drag  force  on  the  ith 
train  is  described  by  g^(x^)  an^  the  applied  force  by  f_^(t). 

As  already  stated,  jerk  (x^)  is  an  important  factor 
affecting  passenger  comfort.  In  order  to  make  jerk  a  variable  in  the 
system  equation  it  is  proposed  to  describe  the  motion  of  the  ith  train 
by  the  third-order  differential  equation 


•  •  • 

m.  x . 
l  l 


+ 


d_ 

dt 


[g.U.)] 


f.(t) 


(2.10) 


It  is  easily  seen  that  2.10  is  obtained  by  differentiating  both  sides 
of  2.9  with  respect  to  time. 
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Since  we  have  assumed  equal  mass  for  all  trains  we  shall 
call  m.=  m  for  i=l, 0  0  ° ,N.  Then  2.10  can  be  written 


g  (x  )  f  (t) 

♦  XI  X 

x„  + -  =  - 

i  m  m 


(2.11) 


A  •  • 

By  defining  x.=  x._,  x„=  x„~  and  x„=  x,„  as  the  state 

l  il  i  i2  l  i3 

variables  we  get  the  following  state  equations: 


Xil(t)  ’  Xi2(t) 


xi2(t)  =  xi3(t) 


(2.12) 


xi3(t)  =  - 


8j(x12)  +  fl(t) 


m 


m 


Now  we  can  define  error  state  variables  as  follows: 


eil(t)  =  xiX(t)  -  z±(t) 


e  9(t)  =  x  (t)  -  v 
i2  i2  s 


(2.13) 


ei3(t)  xi3^t) 


where  z^(t)  is  the  scheduled  position  of  the  ith  train  at  time  t,  vg 
is  the  scheduled  velocity  of  the  ith  train  (the  cruising  velocity) 


and  is  constant. 
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Differentiating  2.13  and  substituting  2.12  yields 


e .,(t)  =  x  (t)  -  z.(t)  =  x  (t)  -  v  =  e  (t)  (2.14) 

11  ll  1  ll  S  1Z 


•  • 


ei2(t)  =  xi2(t)  '  Vs  =  Xi2  =  ei3(t) 


(2.15) 


.  .  ®i^Xi2^  ^^(t) 

e  (t)  =  x  (t)  - - +  - 

i3  i3  mm 


(2.16) 


Taken  together  2.14,  2.15  and  2.16  define  a  new  set  of  state  equations 
describing  the  system. 

As  a  linear  approximation  of  the  drag  force  let 

gi(xi2)  =  d0  +  diei2(t)  (2.17) 

where  d^  is  the  constant  drag  when  the  velocity  of  the  ith  train 
equals  the  scheduled  velocity.  Then 

gi(xi2)  =  =  diei3^  (2.18) 

Now  2,16  becomes 

.  -d  e  ,(t)  +  f . (t) 

e  (t)  =  .  .  - ± -  (2,19) 

m 

•  # 

and  if  we  let  dj^  =  d^/m,  c^  =  f^(t)/u^(t)  and  c|  =  c^/m  where  u^(t) 
is  the  control  variable  for  the  ith  train  we  can  write 

ei3(t)  -  -^ei3(t)  +  cju^t) 


(2,20) 
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Calling  u_^(t)  =  w_^(t)  we  can  write  the  final  state 


equations  for  the  ith  train  as 


eil(t)  =  ei2(t) 


ei2(t)  '  ei3(t) 


(2.21) 


e13(t)  = 


-d^e±3  (t:)  +  c’w.(t) 


Written  in  matrix  form  2.21  becomes 


e’il 

0 

1 

0 

1 

I — 1 

•H 

(D 

1 _ 

0 

ei2 

= 

0 

0 

1 

ei2 

+ 

0 

on  1 

-H 

1 

0 

0 

-di_ 

ei3 

_ci_ 

w±(t) 


(2.22) 


or 


fif  =  As1  +  Bwi(t) 


(2. 23) 


Note  that  w^(t)  is  directly  related  to  the  jerk  forces 


applied  to  the  ith  train.  Specifically, 


H  fi(t) 

w  (t)  =  iL  - — 

1  dt  me! 


f±(t) 


(2.24) 


For  the  sake  of  comparison,  if  the  system  of  equations 


developed  by  Levine  and  Athans  is  written  in  a  form  similar  to  that  of 


_ 
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2.22  we  obtain 


• 

e  •  “1 

0  1 

e . , 

0 

ll 

il 

= 

+ 

• 

e .  _ 

cfl 

1 

O 

e .  _ 

b  . 

i2 

1 

i2 

i 

u.(t) 

1 


where 


a. 

l 


gi(xi2 


'12 


)  - 

Tt) 


fi(t>  -  do 

u±(t) 


(2.25) 


2.4  Cost  Functionals 

Thus  far,  we  have  not  indicated  in  what  manner  we  will 
derive  a  regulator  for  the  system  of  trains.  From  the  discussion  of 
sections  2.1  and  2.2  it  is  evident  that  we  wish  to  minimize  certain 
deviations  from  desired  responses  while  constraining  acceleration,  jerk 
and  headway.  Since  it  is  our  intention  to  use  available  techniques  in 
optimal  control  theory,  it  is  necessary  to  express  our  objectives  in  the 
form  of  a  performance  index  or  cost  functional.  In  this  section  we 
shall  develop  a  general  cost  functional  which  will  express  our  object¬ 
ives  in  a  mathematical  form. 

It  was  mentioned  earlier  that  there  are  two  ways  of 
regulating  a  variable  such  as  say  the  error  in  position  with  respect 

to  scheduled  position,  e#1.  If  we  minimize  e...  for  all  trains  in  the 
r  ll  ll 

system,  each  train  would  need  only  an  on-board  reference  giving  the 
current  scheduled  position  at  all  times.  If,  on  the  other  hand,  we 

minimize  (e..-  e,,.  ,)  for  all  trains  then  each  train  requires  a 

ll  i+l,l 
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knowledge  of  the  position  error  of  all  other  trains  at  all  times. 

Since  the  second  type  of  minimization  would  likely  result  in  a  more 
stable  system  we  shall  concentrate  on  difference  terms  when  setting 
down  cost  functionals. 

The  most  general  cost  functional  considered  is  given  by 


N  2  N_1  2  N"1  2 

Ci-liW±  +i=ii(eil~ei+1»l)  +i^®i(ei2"ei+l,2) 


N 

+  h[( 

i=l 


'il 


N 

+  I  s  !  e 


^  i  i2 


]  dt 


(2,26) 


i=l 


where  N  is  the  total  number  of  trains  and  r, ,  q  .  q! ,  s „  ,  s'  are  non- 

i  i  i  i  i 

negative  weighting  factors. 

The  following  comments  on  2.26  will  explain  the  signi¬ 
ficance  of  the  various  terms  in  the  cost  functional, 

1.  All  terms  in  the  functional  appear  as  quadratics,  so  that 
large  deviations  are  penalized  more  heavily  than  are  small 
deviations.  This  is  a  reasonable  requirement  since  we  wish 
the  system  to  operate  as  close  to  the  schedule  as  possible. 

2.  The  first  term  penalizes  the  system  for  excessive  jerk 
forces  and  hence  reflects  the  passenger  comfort  criteria 
established  in  section  2.2. 

3.  The  second  term  penalizes  the  system  for  differences  in 
position  error  between  successive  trains.  This  has  the 
effect  of  penalizing  the  system  for  deviations  from  the 


scheduled  headway. 
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4. 


5. 


6 . 


7. 


The  third  term  penalizes  the  system  for  differences  in  veloc¬ 
ity  error  between  successive  trains.  At  this  time  it  would 
seem  that  the  addition  of  the  third  term  will  have  little  effect 
on  the  behaviour  of  the  system.  That  is,  if  two  successive 
trains  are  running  with  no  difference  in  position  error  between 
them  theii  their  separation  is  constant.  Hence  they  must  be 
running  at  the  same  speed  so  that  the  difference  in  velocity 
error  between  them  would  be  zero.  If  only  the  third  term  was 
included,  however,  it  would  be  expected  that  the  intervehicle 
spacing  would  not  be  equal  for  all  pairs  of  trains  but  the 
spacing  should  remain  constant. 

The  second  and  third  terms  do  not  guarantee  that  either  the 

position  or  velocity  errors  of  the  individual  trains  will  be 

2  2 

zero.  The  terms  e . and  e..*  L  can  be  large  as  long  as 

il  i  +  l  1 

9 


'(e^-e,  M  .)  remains  close  to  zero, 
il  i+!,l 

To  penalize  the  system  for  deviations  from  the  schedule, 
the  fourth  and  fifth  terms  in  2.24  are  inserted. 

The  upper  limit  on  the  integral  is  specified  as  00  in  order  to 
ensure  that  the  cost  stays  near  zero  after  an  initial  tran¬ 
sient  interval  and  to  avoid  the  specification  of  a  large  term¬ 
inal  time.  Even  if  a  large  terminal  time  were  chosen  we  could 
not  guarantee  that  the  cost  would  remain  near  zero  for  time 
greater  than  the  specified  terminal  time. 

The  cost  functional  is  symmetrical  in  that  positive  and  nega¬ 
tive  jerk  and  deviations  from  position  and  velocity  in  either 


direction  are  penalized  equally. 
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8.  The  constraints  mentioned  in  section  2.2  are  implied  only  in 
as  much  as  an  increase  in  weighting  of  a  particular  term  in 
2.26  has  the  effect  of  constraining  it. 

9.  For  the  second-order  system  described  by  equation  2.25,  the 
general  cost  functional  used  is  very  similar  to  that  of  2.26 
except  the  control  u.  is  substituted  for  w..  That  is. 


JB  ■ 


„  N  ,  N-l  ,  N-l 

l0  1 1  \  + 1  i1(e1rei+i,i)  +  l  si(e12-ei+i>2) 

i=l  i=l  i=l 


+  ^  qieil2  +  ^  Siei22  ^dt  <2'27) 

i=l  i=l 


10.  The  above  cost  functionals  do  not  in  any  way  ensure  good  dyna¬ 
mic  behaviour  of  the  regulator  under  all  possible  conditions 
and  each  system  must  be  simulated  on  a  computer  in  order  to 
check  the  effect  of  the  various  disturbances  on  the  response 
of  the  system. 

In  the  following  chapter  several  different  cost  functionals 
of  the  above  forms  are  used  to  specify  optimal  regulators  for  both  the 
second  and  third-order  representations  of  the  train  system.  In  doing 
so  we  provide  a  wider  basis  of  comparison  between  the  two  representations. 
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CHAPTER  3 


DEVELOPMENT  OF  THE  OPTIMAL  CONTROLLER 


3.1  Method  of  Solution 

Levine  and  Athans  have  shown ^  ^  that  the  problem 
outlined  in  the  preceding  chapter  can  be  solved  by  the  method  that  has 
been  developed  for  the  optimal  control  of  a  linear  system  with  a  quad¬ 
ratic  performance  criterion.  This  problem  is  often  referred  to  as  the 
linear  regulator  problem. 

Although  the  details  of  this  method  are  available  in 
the  literature  and  textbooks  (for  example  see  [7]),  a  summary  of  the 
relevant  results  is  presented  here  for  convenience. 

Consider  a  linear  time-invariant  system  described  by 
the  vector  matrix  differential  equation 

x(t)=Ax(t)+Bu(t)  (3c 1) 

where  x(t)  is  the  n  element  column  state  vector,  u(t)  the  m  element 
column  control  vector  and  A  and  B  are  nxn  and  nxm  constant  matrices 
respectively.  The  system  is  assumed  to  be  controllable  (see  Appendix  1 
for  definition  of  controllability) . 

The  objective  is  to  obtain  a  control  u*(t)  which  will 
take  the  state  x(t)  from  an  initial  state  x(t^)  to  a  final  state  x(T) 
while  minimizing  the  functional 

j  =  JL  /  [<x( t)  ,Qx(t)  >  +  <u(t)  ,Ru(t)  >]  dt  (3,2) 

2  0 

where  Q  is  a  positive  semi-definite  constant  matrix  and  R  is  a  positive 
definite  constant  matrix.  The  control  u(t)  is  unconstrained. 
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It  has  been  shown  that  the  optimal  control  exists,  is 
unique  and  is  given  by 


u*(t)  =  -R  1B'Kx(t) 


(3.3) 


where  the  superscript  prime  denotes  transpose  and  K  is  the  constant 
nxn  positive  definite  matrix  which  is  the  solution  of  the  nonlinear 
matrix  algebraic  equation 


-£A  -  A'K  +  KBR-1B'K  -  Q  =  0 


(3.4) 


The  optimal  trajectory  is  the  solution  of  the  linear  time-invariant 
homogeneous  system 


x(t)  =  Gx(t)  ,  x(0)  given  (3.5) 

where  G  is  defined  by 

G  =  A  -  BR_1B'K 

The  minimum  cost  J*  is  given  by 

J*[x(t)  ]  =  <x(t)  ,Kx(t)  >  (3,6) 

.  [8] 

Concerning  the  matrix  it  has  been  shown  that  the 
assumptions  of  controllability  and  the  terminal  cost  being  zero  imply 
that 

K  =  lim  K(t)  (3.7) 

T-*» 
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where  the  limit  exists  and  is  unique  and  K(t)  is  the  solution  to  the 
Riccati  equation 

K(t)  =  -K(t) A  -  A’K(t)  +  K(t)BR”1B'K(t)  -  Q  (3.8) 

with  boundary  condition  K(T)  =  0  where  T  is  the  final  time. 

Athans  and  Falb  ^  give  the  following  interpretation  to 
3.8.  If  the  time  T  is  thought  of  as  the  "starting  time"  and  K(T)  as 

A 

the  "initial  condition",  then  the  matrix  K  can  be  thought  of  as  the 
"steady-state"  solution  of  the  Riccati  equation  as  the  time  decreases. 

As  T-*»,  the  "transient"  due  to  the  "initial  condition"  dies  down  and 

/\ 

the  solution  approaches  the  value  K  for  t<<T. 

Equation  3.8  can  hence  be  solved  by  a  numerical  method 
by  iterating  backwards  in  time  starting  at  some  large  T.  The  problem 
of  deciding  when  the  steady  state  value  has  been  reached  (i.e.  when 
to  stop  iterating)  is  usually  solved  by  deciding  upon  some  e  such  that 
when  the  results  of  two  successive  iterations  differ  by  e  for  the  first 
time  the  procedure  is  Stopped.  Deciding  on  an  e  can  be  done  by  observ¬ 
ing  the  results  of  the  iterative  scheme  and  in  practice  this  should 
present  no  difficulties. 

3.2  Optimal  Gontrol  For  the  Train  System 

The  method  described  in  the  preceding  section  will  now 
be  applied  to  the  system  of  trains.  We  shall  Carry  out  the  investi¬ 
gation  for  a  string  of  three  trains  so  as  to  keep  the  number  of  equa¬ 
tions  to  a  reasonable  level.  This  choice  will  also  enable  us  to  com¬ 
pare  or  check  our  results  with  those  obtained  by  Levine  and  Athans 


. 
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since  they  also  consider  a  string  of  three  vehicles. 

For  three  vehicles  2.22  becomes 
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This  can  be  written 

e(t)  =  Ag  (t)  +  Bw(t)  (3.10) 

where  A  and  B  are  the  9x9  and  9x3  matrices  of  3.9  respectively.  The 
partitioning  of  3.9  makes  it  easier  to  see  the  coefficients  associated 
with  each  train  in  the  system.  Note  that  3.10  is  of  the  same  form  as 
3.1.  The  cost  functional  in  2.26  can  be  written  as 


1  r  00 

J  =  —  /  [<e(t) ,Qe(t) >  +  <w(t) ,Rw(t) >] dt 

2  0 


(3.11) 
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where 
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0 
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0  0  0  0 
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0  0  0  0 


0  -q2  0  0 

0  0  — s  2  0 

0  0  0  0 


(3.12) 


0  <l2+q3  0  0 

0  0  s2+s3  0 
0  0  0  0 


and 


R  = 


(3.13) 


As  can  be  seen,  both  Q  and  R  are  symmetric.  Observing 
that  the  second,  third,  fourth,  and  fifth  summations  of  2.26  are  >  0  we 
can  say  that  <e(t) ,Qe(t) >>0  hence  Q  is  positive  semi-definite.  Also  it 
can  be  seen  that  R  is  positive  definite. 

Lastly,  we  must  verify  that  the  system  of  3.9  is  completely 
controllable.  We  can  do  this  by  evaluating  S  of  A1.3  (see  Appendix  1) 
for  the  entire  system  or  for  each  train  in  the  siring  (i.e.  show  that 
2.22  is  controllable).  Evaluating  S  for  2.22  gives 


# 
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0 


S 


0  1 


1 


-1 


(3.14) 


It  can  be  verified  by  observation  that  the  3  column  vectors  of  3.14  are 

linearly  independent.  This  also  can  be  verified  by  observing  that 

det  S  =  -1  which  implies  S  is  of  rank  3. 

Hence  we  can  proceed  to  solve  the  Riccati  equation  for 

the  system  of  3.9  obtained  by  using  the  cost  functional  of  3.11. 

Due  to  time  limitations  it  was  decided  to  set 
d.  f  f.(t) 

d!  =  —  =  c!  =  -yi — - —  =  1  for  all  i  and  to  concentrate  on  varying  the 
1  m  1  u_^(t)m 

weighting  factors  within  the  cost  functional  in  order  to  evaluate  per¬ 
formance  of  the  system.  This  is  a  somewhat  arbitrary  decision  but  if 
we  work  on  a  unit  mass  basis  and  take  f^(t)  =  u_^(t)  and  assume  that  the 
drag  force  on  each  train  is  d^  4-  e^CO  then  we  have  d^  =  c|  s  1  which 
will  at  least  let  us  evaluate  regulator  performance. 

With  d!  and  c!  set  equal  to  1  we  can  substitute  into  the 
ii 

Riccati  equation  3.8  for  some  values  of  q4,  ,  r,  ,  s„,  q!  and  s!.  In  choos- 

1  1  1  1  1 

ing  weighting  factors  we  are  really  concerned  with  relative  weighting, 

hence  the  1/2  in  front  of  the  cost  functional  can  be  absorbed  into  the 

weighting  factors  and  r^  set  equal  to  1  for  all  i.  Also,  since  we  wish 

to  penalize  the  system  equally  for  the  misbehaviour  of  any  train,  we 

can  set  q„  =  q ,  q !  =  q ' ,  s .  =  s ,  s !  =  s  for  all  i. 
l  l  i  l 

The  second-order  system  of  2.25  can  be  treated  in  a 
similar  manner  and  a.  and  b.  set  equal  to  1.  We  are  now  in  a  position 

ii 

to  solve  the  Ricatti  equation  for  both  the  second  and  third-order  systems 
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employing  the  following  cost  functionals: 


3  2 

CO  _  O 


h =  K  (t)  +J1{eu(t)-ei+i,i(t)}  ldt 


(3.15) 


3  2 

Jn  =  /  “  [  l  w,2(t)  +5  l  {e  (t)-e  (t) }2 ] dt 

2  0  i=i  .=1  iJ-  1+1»x 


(3.16) 


J„  = 


/ 0°°  [  l  wi2(t)  +  10  l  {eil(t)"ei+l,l(t)} 


]  dt 


i=l 


i=l 


(3.17) 


J4  ’ 


3  .  3 

00  .  r  2 


J  [  l  W  (t)  +  20  l  ^eil(t)_ei+]_  jCO)  ldt 
1=1  1=1 


(3.18) 


JV  = 


3  9  2 

00  -  r>  2 


fo  [  l  w±  (t)  +  5  l  te11(t)-e1+1>1(t)>- 

i=l  i=l 


+  5  l  ^ei2('t')-ei+l  2^^  ■'dt 
i=l 


(3.19) 


J,  =  /  [  l  w  (t)  +  10  l  {e  (t)-e  .  ,  (t)}‘ 

o  o  1  • il  i+l,l 


+  10  l  ^ei2^t^“ei+l  2^^  ■'dt 


i=l 


(3.20) 


j,  - 1:  i 


w42(t)  +  10  l  ^•ei2(-t')“ei+l  2^t')^2^dt 


i=l 


i=l 


(3.21) 


3  0  3  0  2 

°°-r  2  y  x  ,  ,  _  r  2 


Jg  =  Jn  [  l  w_j  (t)  +  10  l  e_.  9  (t)  +  10  l  {en.  n  (t)-e_f+1  -,(t)}  ]dt  (3.22) 


°  i=l  1 


i=l 


'i2 


i=l 


il  i+1,1 


= 


r  t  l  w.2(t)  +10  l  eiJ_2(t)  +10  l  {eil(t)-e1+1#1 
i=l  i=l  i=l 


(t)}  ]dt  (3.23) 


In  the  case  of  the  second-order  system  u_^(t)  is  substituted  for  w_^(t)  in 


the  above. 


1 ' 
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For  each  of  the  cost  functionals,  3.8  can  be  evaluated 
and  put  in  the  form 

K(t)  =  C(t)  (3.24) 

where  C(t)  is  the  right  hand  side  of  3.8.  Equation  3.24  represents 
nxn  simultaneous  nonlinear  first-order  differential  equations  but  since 
K  is  symmetric,  we  need  only  solve  ^  simultaneous  equations. 

Solution  to  2.34  was  obtained  with  the  help  of  an  IBM  360/67  digital 
computer  using  subroutine  RKGS  (fourth  order  Gill's  method)  in  the 
"backwards  mode".  Results  for  the  various  functionals  for  second  and 
third-order  systems  are  shown  in  appendix  2. 

For  the  sake  of  discussing  the  form  of  these  results,  the 
solution  for  the  third-order  system  with  cost  is  shown  below. 
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The  above  matrix  is  partitioned  for  the  purpose  of  the  following  analysis 
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and  can  be  written  as 
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where 


K  = 


K.  .= 


A 
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5*1 
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(3.26) 
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A 

k.  . 
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.j+1 

ki,j+2 

ki+l,j  ki+l,j+l  ki+l, j+2 

AAA 

ki+2,j  ki+2,j+l  ki+2,j+2 


(3.27) 


From  observation  of  3.25  there  are  three  types  of  symmetry  present: 

1.  symmetry  about  the  principal  diagonal 

2.  symmetry  about  the  other  diagonal 

A 

3.  symmetry  of  each  K. . 

[4] 

As  suggested  by  Levine  and  Athans  the  K  matrix  can  be  thought  of  as 
a  "coupling"  matrix  (see  equation  3.3).  Using  this  notion,  the  K^'s 
can  be  interpreted  as  follows : 


K  ^  represents  the  self  coupling  of  train  1 

A 

represents  the  self  coupling  of  train  2 
K 77  represents  the  self  coupling  of  train  3 


*14 

and 

*41 

couple 

the 

1st 

and 

2nd 

trains 

A 

*17 

and 

/A 

*41 

couple 

the 

1st 

and 

3rd 

trains 

A 

*47 

and 

*74 

couple 

the 

2nd 

and 

3rd 

trains 

With  this  interpretation  in  mind  the  three  symmetries  of  3.26  can  be 


explained  as  follows: 
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1.  K  is  symmetrical  (about  the  principal  diagonal)  since  the 
effect  of  say  the  1st  train  on  the  2nd  train  (K^)  is  the 
same  as  the  effect  of  the  2nd  train  on  the  1st  (K^) .  This 
is  reasonable  and  what  we  would  expect. 

2.  K  is  symmetrical  about  the  other  diagonal  since  the  system  is 

A  A 

symmetrical  about  the  middle  train.  For  example  K^= 
because  the  coupling  between  the  1st  and  2nd  trains  is  the 
same  as  that  between  the  2nd  and  3rd  trains.  Again  this  is 
a  reasonable  outcome  and  reflects  the  symmetry  of  the  cost 
functional. 

A 

3.  The  self-coupling  matrices  are  symmetrical  because  the  K 

A  A 

matrix  is  symmetrical.  Further  K  =  K  again  implies  sym- 

n  /y 

metry  about  the  middle  train.  Internal  symmetry  of  the  off- 
diagonal  submatrices  can  be  explained  by  the  following  example. 
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A  A 


17 

CO 

I — 1 

k19 

27 

k28 
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37 

k38 

k39 

(3.28) 


A 

Here  k^g  represents  the  effect  of  the  velocity  of  the  3rd  train 

A 

on  the  position  of  the  1st  train  and  k^y  represents  the  effect 
of  the  position  of  the  3rd  train  on  the  velocity  of  the  1st 

A  A 

train.  The  fact  that  k  Q  =  k..,  indicates  an  internal  symmetry 

lo  Z  / 

between  the  system  state  variables.  Similar  arguments  explain 

A  A  A  A 

why  kig*  kg7  and  k2g=  k3g. 

The  above  discussion  simplifies  the  solution  of  the  Riccati 
equation  for  other  systems  of  this  type  since  the  number  of  equations 


which  need  to  be  solved  is  reduced. 


' 
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3.3  Implementation  of  the  Results 

From  equations  3.3  and  3.5  we  can  construct  the  closed 
loop  optimal  control  system  as  shown  below. 


FIGURE  3-1  THE  CLOSED  LOOP  OPTIMAL  REGULATOR 


From  3.3 


w*(t)  =  -R  1B,Ke(t) 


(3.29) 


Substituting  for  R  and  B  we  obtain 

w*(t)  =  - 

Note  that  due  to  the  nature  of  the  B  matrix  only  the  3rd,  6th  and  9th 
rows  of  K  determine  the  optimal  control  w*(t).  In  a  similar  manner  we 
obtain  the  optimal  control  for  the  second-order  system  as 


^31ell(t)  +  k32e12^  +  •••  +  k39e33^t^ 


k61ell(t)  + 


k91ell(t)  + 


+  k69e33(t) 


+  k99e33^fc^ 


(3.30) 


■ 


■fir 
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u*(t) 


A  A 

k21ell^  +  k22e12^  + 

+  k26e32^t') 

k41ell(t)  +  •  •  • 

+  k46e32(t) 

k61ell(t)  +  •  •  •  +  k66e32 


(3.31) 


In  the  following  chapter  the  optimal  system  is  simulated 
using  both  the  second  and  third-order  representations  and  a  comparison 
of  the  behaviour  under  various  disturbances  is  made. 
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CHAPTER  4  ANALOGUE  COMPUTER  SIMULATION  OF  THE 

TRAIN  SYSTEM 


4.1  Introduction 

Using  the  solutions  to  the  Riccati  equations  corresponding 
to  the  various  cost  functionals  obtained  in  the  previous  chapter  and 
using  the  closed  loop  configuration  shown  in  figure  3-1  both  the  second 
and  third-order  models  of  the  train  system  were  simulated  on  a  PACE  231R 
analogue  computer.  The  analogue  computer  diagrams  for  the  second  and 
third-order  models  are  shown  in  figures  4-1  and  4-2  respectively.  In 
this  chapter,  comparisons  in  the  performance  of  the  optimal  control 
systems  for  the  second-order  and  third-order  models  using  the  different 
cost  functionals  are  presented. 

In  order  to  assess  the  usefulness  of  the  two  models  they 
were  both  subjected  to  the  same  series  of  tests  which  are  outlined  below: 

1.  Train  //I  is  initially  going  too  fast.  This  condition  is 
simulated  by  setting  positive  initial  conditions  on  the 
velocity  and  position  errors  for  the  first  train.  This 

test  was  performed  on  both  the  second  and  third-order  systems 

for  the  cost  functionals  J, ,  J„,  J,,  Jr  ,  J-  and  J„  defined 

1  3  4  6  7  9 

in  chapter  3.  In  the  work  of  Levine  and  Athans^^]  only 
the  simulation  of  the  second-order  train  system  is  considered 
and  that  for  only  one  cost  functional  (J^) . 

2.  A  small  step  change  in  the  velocity  error  of  train  #2 .  This 
is  the  most  interesting  test  in  so  far  as  comparing  the 
second  and  third-order  systems  is  concerned  in  that  the 
relationship  between  the  respective  jerk  and  acceleration 


'll 
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Figure  4-1  ANALOGUE  SIMULATION  OF  SECOND-ORDER  SYSTEM 


Figure  4-2  ANALOGUE  SIMULATION  OF  TRAIN  #1  FOR 
THIRD-ORDER  SYSTEM 
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forces  is  most  apparent. 

3.  Catastrophic  failure.  A  negative  ramp  change  in  velocity 
error  of  train  //I  was  superimposed  on  the  system  to  simulate 
the  first  train  decelerating  at  a  high  and  uncontrolled  rate. 
Before  discussing  the  results  of  the  above  tests  it  may 
be  useful  to  indicate  the  relationships  between  the  second  and  third- 
order  control  variables  and  the  acceleration  and  jerk  forces. 

For  the  second-order  system,  from  2.25 


u±(t) 


f . 

l 


( t )  —  d, 


0 


b. 

i 


(4.1) 


where  d^  is  the  drag  force  on  the  ith  train  when  there  is  no  velocity  error 
and  b_^  is  a  constant  of  proportionality. 

By  differentiating  4.1  with  respect  to  time  we  get 


u.(t) 


f  ±  (t ) 


b. 

l 


(4.2) 


which  relates  the  control  variable  to  jerk  forces  for  the  second-order 
system. 

For  the  third-order  system,  from  2.24 


w±(t) 


fi(t) 


c. 

i 


(4.3) 
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and  integrating  the  above  yields 

,  00  f  .  (t) 

/  w.(T)dx=_i -  (4.4) 

0  1  c. 

1 

which  relates  the  accelerating  force  to  the  control  variable  for  the 
third-order  system. 

4.2  Results  of  Test  1 

In  the  following,  we  will  compare  the  overall  behaviour 

of  the  second  and  third-order  control  systems  as  well  as  the  effect  of 

the  different  cost  functionals  on  the  two  systems. 

The  ratio  of  e^^(O) /e^2 (0)  was  chosen  arbitrarily  to  be 

10/3  for  all  the  following  tests.  Referring  to  section  3.2  it  can  be 

seen  that  cost  functionals  J1 ,  J  and  J,  are  of  the  same  form  but  have 

2  2 

different  weighting  on  the  term  J  (e  -  e  -i)  which  penalizes  for 

i  1  i  '  1 » I 

error  in  headway  from  schedule.  Second-order  system  response  to  test  1 
using  these  cost  functionals  is  shown  in  figures  4-3(J^),  4-4(J2) ,  and 
4-7(J^).  The  corresponding  third-order  system  response  is  shown  in 
figures  4-5  and  4-6(J2)  and  4-10(J^,  and  J^) .  It  should  first  be 
noted  that  the  second-order  system  response  for  cost  agrees  with  that 
obtained  by  Levine  and  Athans  which  is  some  assurance  that  our  simu¬ 
lation  is  correct. 

First,  it  may  be  useful  to  note  the  similarities  in  res¬ 
ponse  of  the  second  and  third-order  systems  using  these  cost  functionals. 
1.  As  predicted  in  section  2.4,  the  cost  functionals  J^,  and 
do  not  guarantee  a  zero  steady  state  position  error  for 
any  of  the  three  trains.  However,  the  steady  state  value  of 


1 


: 
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position  error  for  all  trains  is  equal,  constant  and  less  than 
the  value  at  time  t=0.  The  actual  value  at  steady  state  varies 
depending  on  the  cost  functional  used  (relative  magnitudes  are 
indicated  in  the  figures) . 

2.  For  both  systems,  with  increased  weighting  on  the  term 


2 


2 

e^+^  in  the  cost  functional,  the  attempted  correction 


in  position  error  of  a  given  train  takes  place  in  less  time  and 
is  accompanied  by  a  greater  amount  of  overshoot  in  the  state 


variables  e._  and  e._.  Excessive  overshoot  of  any  state  variable 
ll  12 


is  undesirable  from  the  standpoint  of  passenger  comfort.  Approx¬ 
imate  transient  time  intervals  for  each  system  are  indicated  on 
the  figures. 


Next,  it  is  possible  to  point  out  the  differences  in  behav¬ 


iour  of  the  second  and  third-order  systems  for  equivalent  cost  functionals. 
1.  Generally,  the  transient  intervals  for  the  second-order  systems 
are  less  than  those  for  the  corresponding  third-order  systems. 


For  instance,  in  the  case  of  cost  J^,  the  second-order  system 


(figure  4-4)  exhibits  about  5  seconds  of  transient  behaviour 
while  the  third-order  system  (figure  4-5)  exhibits  a  transient 
time  of  about  9  seconds. 

2.  The  damping  of  the  second-order  system  for  a  given  cost  func¬ 
tional  appears  to  be  greater  than  that  for  the  corresponding 
third-order  system.  For  example,  considering  the  systems  with 
cost  (figures  4-4  and  4-5)  ,  the  state  e^(t)  clearly  shows 
less  overshoot  in  the  case  of  the  second-order  system  than  it 


does  in  the  third-order  system. 


. 


V 
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3.  The  value  of  steady  state  position  error  for  the  trains,  given 
any  one  of  the  three  cost  functionals  considered,  is  less  for 
the  third-order  system  than  for  the  second-order  system,, 

Let  us  now  consider  the  systems  corresponding  to  the  cost 
functionals  and  J^.  Referring  to  figures  4-8(J^)  and  4-9 (J^)  for 
the  second-order  system  and  to  figure  4-10  for  the  third-order  system, 
the  following  comments  can  be  made. 

1.  Cost  functional  is  identical  to  except  for  the  added 

t  2 

term  10  l  (e^-  ei+l  2)  penalizes  the  system  for 

i=l  * 

differences  in  velocity  errors  between  successive  trains. 

As  was  predicted  in  section  2.4,  the  addition  of  this  term  does 
not  appreciably  affect  the  behaviour  of  either  the  second  or 
third-order  systems  compared  to  that  using  cost  Jy  Comparing 
the  second  and  third-order  systems  using  cost  ,  one  can 
make  the  same  observations  as  were  made  concerning  the  systems 
using  cost  J^,  and  J^. 

2.  Cost  functional  penalizes  the  system  only  for  differences 

in  velocity  errors  between  successive  trains  and  as  can  be 
seen  in  figure  4-9  (second-order  system)  the  steady  state 
position  errors  are  different  for  each  train  in  the  string. 

This  is  to  be  expected  since  no  penalty  is  imposed  for  differ¬ 
ences  in  position  errors  between  successive  trains  as  was  the 

case  for  cost  functionals  J, ,  J_,  J,  ,  and  J, .  Hence  this  con- 

1  3  4’  6 

troller  does  not  regulate  for  headway  and  cannot  be  considered 
satisfactory  in  meeting  our  objective.  Although  only  the 
state  e^(t)  s^own  f°r  ^7  f°r  the  third-order  system  in 
figure  4-10,  system  behaviour  is  very  similar  to  that  shown 
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in  figure  4-9  for  the  second-order  system. 

Lastly,  let  us  consider  the  systems  corresponding  to  cost 

3  2 

Jq  which  differs  from  cost  J-  in  that  the  term  10  £  e._  is  added  inside 

J  i=l  11 

the  integral.  This  has  a  marked  effect  on  the  system  behaviour  to  test  1 
in  that  steady  state  position  errors  for  all  trains  are  zero.  Referring 
to  figures  4-11  (second-order)  and  4-12  (third-order)  the  following  can 
be  stated. 

1.  Since  steady  state  position  errors  are  zero  for  all  trains, 
this  controller  regulates  for  headway  and  scheduled  position. 

Also,  for  both  the  second  and  third-order  systems  the  initial 
condition  of  train  #1  going  too  fast  does  not  disturb  the  fol¬ 
lowing  trains  to  a  great  extent  during  the  correction  interval. 

This  is  of  course  due  to  the  fact  that  all  position  errors 
remain  at  zero  for  the  following  trains  once  the  transient 
interval  has  passed. 

2.  Comparing  the  two  systems,  one  can  easily  see  that  the  tran¬ 
sient  period  is  greater  for  the  third-order  system.  Both 
systems  display  some  overshoot  which  could  be  eliminated  by 
reducing  the  weighting  on  both  error  terms  in  the  cost  func¬ 
tional  Jg .  Generally,  this  cost  functional  results  in  the 
controller  best  meeting  our  objectives  for  both  systems.  As 
a  result  of  test  1,  we  cannot  say  which  of  the  two  systems  is 
to  be  preferred. 

4.3  Results  of  Test  2 

To  simulate  a  negative  step  change  in  velocity  error  (e^) 
of  the  first  train,  a  negative  step  voltage  change  in  e ^  was  automatically 
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Figure  4-3.  STATE  RESPONSE  TO  TEST  1  OF  2ND-ORDER  SYSTEM  WITH  COST  J 
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Figure  4-4  STATE  RESPONSE  TO  TEST  1  OF  2ND-ORDER  SYSTEM  WITH  COST  J3 
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Figure  4-10  POSITION  ERROR  OF  TRAIN  1  VS.  TIME  FOR  THIRD-ORDER 
SYSTEM  WITH  VARIOUS  COST  FUNCTIONALS. 
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switched  into  the  circuit  after  3  seconds  of  unperturbed  operation. 

Cost  functional  was  used  in  this  test.  Figures  4-13a  and  4-13b  show 
the  second-order  system  control  variable  and  its  derivative  as  a  function 
of  time  (for  the  first  train).  The  step  change  in  e^tt)  occurs  at  t=3 
seconds  at  which  time  the  control  u^(t)  undergoes  a  step  change  (figure 
4-13a)  and  then  decays  to  zero  again  in  about  4  seconds.  The  step  change 
in  u^(t)  at  t=3  seconds  results  in  a  step  change  in  accelerating  force 
(see  equation  4,1)  and  hence  an  infinite  jerk  force.  This  is  verified 
in  figure  4-13b  though  the  pen  recorder  is  unable  to  exactly  follow  the 
response.  The  noise  on  this  plot  is  due  to  the  differentiator  used  and 
was  reduced  somewhat  by  use  of  a  low-pass  filter. 

The  third-order  system  was  subjected  to  the  same  step 
change  in  e ^  as  the  second-order  system  and  the  results  are  shown  in 
figures  4-13c  and  4-13d.  The  third-order  system  control  for  the  first 
train,  w^(t),  behaves  in  a  similar  manner  as  does  u^(t)  in  the  case  of 
the  second-order  system.  However,  since  from  4.3,  w^(t)  is  directly 
proportional  to  jerk  force,  no  infinite  jerk  is  experienced  as  was  the 
case  for  the  second-order  system.  Comparing  figures  4-13d  and  4-13a 
the  accelerating  forces  of  the  two  systems  can  be  compared.  Although 
the  accelerating  force  is  applied  to  the  third-order  system  for  a  greater 
length  of  time  than  in  the  case  of  the  second-order  system,  the  appli¬ 
cation  of  the  force  is  made  in  a  smoother  manner,  eliminating  infinite 
jerk. 

The  differences  in  the  two  systems  indicated  above  will 
be  evident  for  any  sudden  change  in  one  of  the  states  since  the  accel¬ 
erating  force  for  the  third-order  system  is  the  integral  of  the  control 
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variable  while  for  the  second-order  system  the  accelerating  force  is 
directly  proportional  to  the  control  variable, 

4.4  Results  of  Test  3 

In  the  previous  two  tests  no  attempt  was  made  to  assign 
values  to  the  position,  velocity,  acceleration  and  jerk  curves  although 
this  could  be  done.  In  the  third  test,  in  which  the  first  train  was 
simulated  to  be  decelerating  at  an  uncontrolled  rate,  we  have  attempted 
to  make  some  quantitative  statements  about  the  behaviour  of  the  systems. 

To  simulate  a  constant  deceleration  of  train  //I  the  feed¬ 
back  paths  for  e^  and  were  opened  and  a  negative  ramp  voltage  sub¬ 
stituted  for  e^  and  its  integral  for  e^.  The  results  °f  the  test 
using  cost  functionals  JQ  and  J„  for  both  the  second  and  third-order 
systems  are  shown  in  figures  4-14,  4-15,  4-16  and  4-17.  By  differen¬ 
tiating  the  second  equation  of  2,13  we  get 

ei2(t)  =  xi2(t)  (4.5) 

and  hence  the  slope  of  the  imposed  e^2(t)  ramp  can  be  assinged  some 

deceleration  value  which  will  be  equal  to  the  deceleration  rate  of  train  //l. 

2 

For  this  test  the  value  of  1.0  g.  (or  32.2  ft. /sec.  )  was  arbitrarily 
given  to  e^2(t)  '  ^nce  this  value  is  established,  the  position  states 
can  be  calibrated  in  feet  and  the  velocity  states  calibrated  in  ft. /sec. 

Under  constant  deceleration,  a  collision  between  trains 
1  and  2  will  occur  when  for  the  first  time 


hr  e2.i.l  =  ho 


,  k>0 


(4.6) 
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Figures  4-18  (second-order)  and  4-20  (third-order)  show  the  relationship 
between  |e^-  e2l|  anc*  f°r  both  cost  functionals  Jg  and  Jg.  From 

4.6,  at  collision,  Je.^-  e2l|  =  ^0*  t^e  ^eac^way»  so  these  figures  can 
be  thought  of  as  plots  of  headway  vs  time  before  collision  between  trains 
1  and  2.  For  a  given  collision  time  we  can  then  determine  the  velocity 
of  impact  given  by 


Vimpact  e22  e12 


(4.7) 


Figures  4-19  and  4-21  are  plots  of  headway  vs  v^mpact  f°r  the  second 
and  third-order  systems  respectively.  The  following  comments  on  test  3 
can  now  be  made. 

1.  For  both  the  second  and  third-order  systems,  cost  resulted 
in  better  maintenance  of  spacing  between  the  1st  and  2nd 
trains  than  did  cost  JQ  (as  evidenced  in  figures  4-18  and 

O 

4-20).  This  is  due  to  the  penalty  that  cost  Jg  imposes  on  the 
system  for  large  velocity  errors.  At  any  time,  t,  the  velo¬ 
city  error  of  train  2  with  cost  JQ  is  less  than  that  with 

O 

cost  J~  which  means  that  the  system  with  cost  JQ  is  unable 
3  o 

to  make  corrections  in  the  position  of  the  2nd  train  as 
rapidly  as  the  system  with  cost  Jg. 

2.  In  comparing  second  and  third-order  system  behaviour,  it 
should  be  noted  that  in  the  case  of  the  second-order  system 
with  both  cost  functionals  Jc  and  J„,  for  large  enough  head- 
way  the  impact  velocity  approaches  a  constant  (figure  4-19). 

In  other  words,  there  is  a  small  delay  on  the  part  of  the  2nd 
train  in  reaching  the  deceleration  rate  of  the  1st  train  and 
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thus  for  large  headways  (corresponding  to  high  cruising 
velocities)  the  delay  causes  an  initial  offset  in  (e^-  e12^ 
after  which  both  trains  decelerate  at  the  same  rate. 

For  the  third-order  system  with  cost  JQ,  the  decel- 

O 

eration  rate  of  the  2nd  train  never  becomes  as  great  as  that 
of  the  1st  train  and  hence  the  impact  velocity  is  proportional 
to  headway  for  large  enough  headway  (figure  4-21) . 

For  the  third-order  system  with  cost  J^»  the  2nd  train 
decelerates  at  a  greater  rate  than  the  1st  train  thereby  over¬ 
coming  the  time  delay  effect.  After  about  13  sec.  (figure  4-17) 
both  trains  decelerate  at  the  same  rate.  As  can  be  seen  in 
figure  4-21  there  will  be  no  collisions  for  headways  greater 
than  300  feet  (for  a  deceleration  rate  of  1.0  g.).  Due  to 
time  limitations,  this  unusual  behaviour  of  the  third-order 
system  was  not  investigated  further. 
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Figure  4-15  RESULTS  OF  TEST  3  FOR  3RD-ORDER  SYSTEM  WITH  COST  J0 
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Figure  4-16  RESULTS  OF  TEST  3  FOR  2ND-ORDER  SYSTEM  WITH  COST  J 
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Figure  4-17  RESULTS  OF  TEST  3  FOR  3RD-ORDER  SYSTEM  WITH  COST  J3 
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Figure  4-18  /  e..  -  e../VS.  TIME  FOR  SECOND-ORDER  SYSTEM  SUBJECT  TO  TEST  3 
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Figure  4-19  HEADWAY  VS.  VELOCITY  OF  IMPACT  FOR  THE  SECOND-ORDER  SYSTEM 
SUBJECT  TO  TEST  3 
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Figure  4-20  /e.,  -  e../ VS.  TIME  FOR  THIRD-ORDER  SYSTEM  SUBJECT 
TO  TEST  3 
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Figure  4-21  HEADWAY  VS.  VELOCITY  OF  IMPACT  FOR  THE  THIRD-ORDER  SYSTEM 
SUBJECT  TO  TEST  3 
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CHAPTER  5  SUMMARY,  CONCLUSIONS  AND  SUGGESTIONS 

FOR  FUTURE  RESEARCH 

5.1  Summary 

Before  making  any  conclusions  it  will  be  useful  to 
briefly  summarize  what  has  been  done  in  this  thesis  in  so  far  as  the 
development  of  optimal  controllers  is  concerned.  First,  the  second- 
order  mathematical  model  developed  by  Levine  and  Athans  was  studied  in 
some  detail  by  simulation  on  the  analogue  computer.  A  greater  number 
of  cost  functionals  was  utilized  than  in  the  previous  study  of  this 
model  thereby  enabling  evaluation  of  the  system  under  a  greater  variety 
of  constraints.  Second,  the  third-order  model  of  the  train  system 
developed  in  this  thesis  was  simulated  and  its  performance  under  a 
variety  of  tests  was  compared  to  that  of  the  second-order  system. 

5.2  Conclusions 

On  the  basis  of  the  results  reported  in  the  preceding 
chapter  the  following  conclusions  may  be  stated. 

1.  The  response  of  the  train  system  to  any  perturbation  depends 

to  a  great  extent  on  the  index  of  performance  or  cost  functional 
chosen.  This  is  of  course  to  be  expected  since  our  only  crite¬ 
rion  determining  the  behaviour  of  the  system  is  the  minimization 
of  the  cost  functional  over  the  time  interval  during  which  the 
control  variable  is  non-zero.  We  can  therefore  exert  a  great 
deal  of  control  over  the  regulator  performance  by  altering  the 
weighting  factors  of  the  terms  within  the  cost  functional. 
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2.  To  achieve  an  optimal  controller  which  will  regulate  for  both 
headway  and  position  schedule  it  is  necessary  to  specify  a  cost 
functional  of  the  form 


negative  weighting  factors. 

This  statement  is  true  for  both  the  second  and  third- 
order  models . 

3.  On  the  basis  of  the  tests  conducted,  the  third-order  model 
appears  to  result  in  a  train  regulator  which  better  meets  the 
requirement  that  excessive  jerk  forces  be  minimized.  More 
research  on  the  reaction  of  the  third-order  system  to  various 
types  of  catastrophic  failure  would  be  needed  before  any  con¬ 
clusions  could  be  made  concerning  this  phase  Qf  the  study.  In 
the  interest  of  safety  and  reliability,  an  overriding  control 
system  that  would  function  to  apply  maximum  stopping  power  to 
all  trains  once  any  train  was  detected  to  be  decelerating  at 
an  abnormally  high  rate  would  probably  be  most  satisfactory. 
Under  this  condition  equation  2.8  would  govern  the  number  of 
collisions . 

4.  It  should  be  clear  from  the  work  of  this  thesis  that  it  is 
possible  to  thoroughly  analyze  a  train  regulating  system  by 
developing  a  mathematical  model  of  the  string  of  trains,  apply¬ 
ing  known  control  theory  to  develop  a  controller,  and  simulating 
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the  resulting  system  on  an  analogue  computer.  Hence  it  is 
unnecessary  to  use  trial  and  error  approaches  which  have 
characterized  the  area  of  train  technology  in  the  past. 

5.3  Suggestions  for  Future  Research 

Possible  future  research  on  the  subject  of  control  of 
high  speed  vehicles  falls  into  two  categories  —  straightforward 
extensions  of  the  present  work  and  analysis  of  completely  different 
aspects  of  high  speed  train  operation. 

Under  the  first  category  it  would  be  profitable  to 
re-examine  the  optimal  controllers  obtained  when  the  constants  d^  and 
c^  in  2.22  are  set  to  values  other  than  unity.  Increasing  values  for 
d^  and  c^  would  have  the  effect  of  reducing  the  mass  and  hence  inertia 
of  the  ith  train. 

As  previously  indicated  we  can  "constrain"  certain 
variables  in  the  system  only  by  altering  the  weighting  terms  in  the 
cost  functionals.  It  would  be  useful  to  determine  an  optimal  regulator 
using  other  theory  of  optimal  control  (such  as  non-linear  programming 
or  Pontryagin’s  Maximum  Principle)  in  which  constraints  such  as 
J w  ( t ) | < 1  can  be  considered.  This  might  require  a  different  formu¬ 
lation  of  the  mathematical  model  (such  as  discreet  differential  equations 
in  the  case  of  non-linear  programming)  but  might  lead  to  a  more  general 
result . 

Under  the  second  category  of  future  research  is  the 

[9] 

problem  of  the  merging  of  two  strings  of  vehicles.  Athans  has 
reported  the  following  approach.  The  physical  problem  is  indicated  in 


Figure  5-1. 
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FIGURE  5-1  THE  MERGING  PROBLEM 


The  problem  can  be  broken  into  two  parts  —  first,  what 
is  the  best  merging  order  and  second,  given  a  merging  order  what  is  the 
best  way  to  control  acceleration  and  velocity?  If  it  is  assumed  that 
merging  begins  far  back  from  the  merge  point  then  for  a  given  merging 
sequence  the  vehicles  can  be  considered  to  be  on  one  track.  The  number 
of  possible  merging  sequences  is  given  by 


=  (q+g) ! 

a !  0  ! 


(5.2) 


So  for  the  yth  possible  merging  sequence,  the  motion  of 
the  ith  vehicle  is  governed  by  the  differential  equation 


•  ♦ 


X  . 

iy 


(t) 


-  -ax.  (t) 
i  iy 


biuiY(t) 


(5.3) 


where  x(t)  is  the  position  of  the  ith  train  and  a^  and  bi  are  the  drag 
and  mass  coefficients  of  equation  2.25. 


y  si"  >  }» rl .  -i  i  -  3  f  !t>  :  rt  i  £  ~t  j  J  :  <1  ft 

_ 
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Error  state  variables  can  be  introduced  as  before  giving 
for  the  y th  merging  sequence  the  system  of  2.25.  The  acceleration  and 
velocity  are  controlled  by  assuming  cost  functionals  indentical  to  those 
of  section  2.4  for  each  sequence.  For  each  merging  order  we  can  com¬ 
pute  the  cost 


Jy  =  j  <ey(t),  Kgy(t)>  (5.4) 

and  chose  y*  such  that  J^<  for  all  Y^Y*j  hence  determining  the  best 
merging  order. 

It  is  thought  that  since  the  third-order  formulation 
yielded  a  better  regulator  from  the  standpoint  of  eliminating  excessive 
jerk,  it  might  be  profitable  to  use  this  formulation  in  place  of  that 
above  and  compare  merging  performance  of  the  two  systems. 

If  the  third-order  system  proved  satisfactory  in  this 
regard,  the  optimal  control  system  could  then  regulate  most  of  the 
necessary  movements  of  an  actual  system  of  trains  including  slowing 
down  for  curves,  stations  etc.  and  re-merging  with  the  main  guideway. 
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APPENDIX  1 


THE  MEANING  OF  STATE-CONTROLLABILITY 


y 


m  inputs 


I 

p  outputs 


1> 


n  states 


V 


a 


A  system  G  defined  by  the  state  equation 

x(t)  =  Ax(t)  +  Bu(t)  (Al.l) 

is  said  to  be  state-controllable  if  knowing  the  matrices  A  and  B  it  is 
possible  to  construct  an  input  u  which  will  take  the  state  x(t^)  to 
x(tf)  =  0  in  some  finite  time  interval  T. 

It  is  shown  in  [7]  that  the  initial  state  can  be  written 
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where 

Yk(t0>T)  =  “  /  f  otk(t0-x)y(x)dT  (A1.4) 

From  A1.3, 

V  =  flx(tQ)  (A1.5) 

which  by  A1.4  stipulates  u(t).  Hence  A1.5  must  be  valid  which  means 
that  S  ^  must  exist  (i.e.  det  S^O) .  This  requires  that  S  (dimensions 
nxnm)  be  of  rank  n  for  our  system  to  exhibit  state  controllability,  [n  is 
the  dimension  of  x(t)  and  m  the  dimension  of  u(t)] 


(A.XA) 

>j  Jw  bi  lbv  ad  Jeum  2.XA  9DnaH  .  (x)y  asXBit/qiXa  A.XA  vd  rloiriw 

j.  .  Lc  'J  :>  sx  a  dln'xfixa  cX  n  -jJeve  cct.  xa  sc  (  iaxn 
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